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Abstract: We compute the complete next-to-leading order (NLO) SUSY-QCD correc- 
tions for the associated production of a pseudoscalar Higgs boson with a bottom quark 
via bottom-gluon fusion at the CERN Large Hadron Collider (LHC) and the Fermilab 
Tevatron. We find that the NLO QCD correction in the MSSM reaches 40% ~ 50% at the 
LHC and 45% ~ 80% at the Tevatron in our chosen parameter space. 
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1. Introduction 

The standard model (SM) ]lj], with the Higgs mechanism manifested by a single complex 
scalar weak isospin doublet that spontaneously break the electrowak gauge symmetry, 
is believed to be an incomplete description of nature. There exists the problem of the 
quadratically divergent contributions to the corrections to the Higgs boson mass. This is 
the so-called naturalness problem of the SM. One of the good methods to solve this problem 
is to make supersymmetric (SUSY) extensions to the SM. The quadratic divergences of the 
Higgs mass can be cancelled by loop diagrams involving the supersymmetric partners of 
the SM particles exactly. The most attractive supersymmetric extension of the SM is the 
minimal supersymmetric standard model (MSSM) |3|]. Any enlargement of the sector 
beyond the single SU(2)l Higgs doublet of SM, with two or more doublets as required in 
supersymmetric theory, necessarily involves new physical particles. In the MSSM, there 
are two Higgs doublets H® and to give masses to down- and up- type fermions. The 
Higgs sector consists of three neutral Higgs bosons, one CP-odd particle (A ), two CP-even 
particles {hP and H°), and a pair of charged Higgs bosons (H^). 

The discovery of the additional heavy Higgs bosons will help us to probe the contents 
of the supersymmetric Higgs sector. Until now all supersymmetric Higgs bosons haven't 
been directly explored yet, only LEP2 group presents the strongest lower mass limits of 
91.0 GeV and 91.9 GeV for the light CP-even and the CP-odd neutral Higgs bosons h° and 
respectively. The CERN Large Hadron Collider (LHC), which is a proton-proton 
collider with ^/S = 14 TeV and a luminosity of 100 fb _1 per year, has been designed 
specifically to continue finding Higgs bosons. In the MSSM theory with a large value of 
tan/3(tan/3 = V2/V1, v\ and V2 are the vacuum expectation values of the two Higgs boson 
fields Hf and H®, respectively), the strength of the A°bb coupling increases greatly. The 
pseudoscalar Higgs boson ^4° could be produced with a substantial rate at the LHC either 
in association with bottom quarks [qq,gg — ► bbA° or gb(b) — > b(b)A°] or through gg — » ^4° 
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when b or b loops dominate, provided that tan/3 is large enough||. Ref.|| presents the 
calculation of the total cross section of the ^4° plus two jets at the lowest order. The total 
cross section for the inclusive production of A has been calculated at NLO in Ref.Q and 
at NNLO in Refs. ||] || , where the authors use the effective lagrangian for the interaction 
of ^4° Higgs boson with the gluons and neglect the contribution of bottom quark loop so 
that their result is only for small and moderate value of tan (5. 

Because the high-pr bottom quark can be tagged with reasonably high efficiency, the 
observation of a bottom quark with high px can reduce the backgrounds of the A Higgs 
boson production. The A Higgs boson production associated with bottom quark can occur 
via tree-level subprocess gb(b) — > b(b)A°[l3\, where the initial bottom quark resides in the 
proton sea. The cross section for the production of Higgs boson h° and a single high px 
bottom quark via subprocess gb(b) — > b(b)h° has been studied in both QCD and SUSY- 
QCD at NLO [14, [D], 16]. However, it is pointed out in Ref.[10] that the calculation of 



bb — ► h° may overestimate the inclusive cross section with the introduction of conventional 
b densities, due to crude approximations inherent in the kinematics, which give rise to 
large bottom quark mass and phase space effect. But in Ref.[ll], it is shown that the 
bottom parton approach is valid by choosing the appropriate factorization scale for the 
process bb — > h° with fJ-f = mh/^ rather than m^. In Ref.jjl^], the author computed 
the NLO contributions to the inclusive cross section pp —* tH~ via the subprocess bg — > 
tH~ and shown the bottom parton approach is valid for this process by choosing the 
factorization scale with about [if ~ m av /3 = (m t + m^-)/6. In our calculation, we choose 
the factorization scale as fj,f = m^/i when we use the bottom parton approximation. 

In this paper, we calculated the cross section for the associated production of the ^4° 
Higgs boson and a single high-p^ bottom quark via bg — ► A°b(hg — > A%) in the MSSM at 
the Tevatron and LHC including the NLO QCD corrections. The structure of this paper 
is as follow: In Sec. 2, we discuss the LO results of the subprocess bg — > A°b. In Sec. 3, 
we present the calculations of the NLO QCD corrections. In Sec. 4, the numerical results 
and conclusions are presented. 



2. The leading order cross section 

Since the cross sections for the subprocess bg — > A°b and its charge-conjugate subprocess 
bg -> A°b in the CP-conserved MSSM are same, we present only the calculation of the sub- 
process b(pi)g{p2) — ► A \kz)b(ki) here (where pi^ and ^4 represent the four-momentum of 
the incoming partons and the outgoing particles, respectively). The subprocess bg — > A°b 
can occur through both s-channel and t-channel as shown in Fig.l(A-B). So we divide the 
tree-level amplitude into two parts and denote it as 

M° = M^ ) + M^ t) , (2.1) 

where and represent the amplitudes arising from the s-channel diagram shown 
in Fig. 1(A) and the t-channel diagram shown in Fig. 1(B) at the tree- level, respectively. 
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(A) 

Figure 1: Leading order Feynman diagrams for the subprocess of bg — » A°b 
The explicit expressions for the amplitudes Mq and M^f 1 can be written as 

M (s) = ^^)^^) iZ . (fc4 ) 75 (^ 1+ ^ 2 ) 7 ^. (pi ) e g (p2) ^ ) 

M? = ^^^ (frh^l _ ^ )75U . (pi)e a (p2) ^ ) (2 . 2) 

where s = (pi +P2) 2 , t = {pi — k^) 2 and u = {p\ — k^) 2 are the usual Mandelstam variables. 
fj, r is the renormalization scale, g s (fi r ) is the running strong coupling strength and T a is 
the SU(3) color matrix. VAbb{l^r)lb is the Yukawa coupling between A Higgs boson and 
bottom quarks. In MSSM, V^bi^r) is given as 

m&(// r ) is the MS running mass of the bottom quark. It is well known that we should 
use the running mass rather than the pole mass when evaluation of the Yukawa coupling, 
because the pole mass Yukawa coupling will yield a huge overestimate of the cross section. 
Throughout our evaluation we neglect the bottom quark mass except in the Yukawa cou- 
plings. This corresponds to the simplified Aivazis-Collins-Olness-Tung(ACOT) scheme [|l9|| . 
In any diagram in which the bottom quark is an initial-state parton, the bottom quark 
mass may be neglected without any loss of accuracy. 

Then the lowest order cross section for the subprocess bg — > A°b in the MSSM is 
obtained by using the following formula: 

3°{6, bg - AH) = ^ f maX di |>T, (2.4) 

where i max = and i m in = m \ ~ The summation is taken over the spins and colors 
of initial and final states, and the bar over the summation recalls averaging over the spins 
and colors of initial partons. 



3. NLO QCD corrections 

The NLO QCD contributions to the subprocess bg — > A% can be separated into two parts: 
the virtual corrections arising from loop diagrams and the real gluon emission corrections. 
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3.1 Virtual Corrections 



The virtual corrections in the MSSM to bg — > A°b consist of self-energy, vertex and box 
diagrams which are shown in Figs. 2-3. Fig. 2 shows the one-loop diagrams of the SM-like 
QCD corrections from quarks and gluons, and Fig.3 presents the one-loop diagrams of 
the SUSY QCD corrections from squarks and gluinos. There exist both ultraviolet (UV) 
and soft /collinear infrared(IR) singularities in the amplitude from the SM-like diagrams in 
Fig.2, and the amplitude part from SUSY QCD diagrams (Fig. 3) only contains UV singu- 
larities. In our calculation, we adopt the 't Hooft-Feynman gauge and all the divergences 
are regularized by using dimensional regularization method in d = 4 — 2e dimensions. 




(Q) 



Figure 2: Virtual one-loop Feynman diagrams for the subprocess of bg — > A°b of the SM-like QCD 
corrections. 

In order to remove the UV divergences, we need to renormalize the wave functions 
of the external fields, the strong coupling and the A — b — b Yukawa coupling. For the 
renormalization of the strong and Yukawa couplings, we employ the modified Minimal 
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Figure 3: Virtual one-loop Feynman diagrams for the subprocess of bg — > A°b of the SUSY QCD 
corrections. 

Subtraction (MS) scheme. The relevant renormalization constants in this work are defined 
same as those in Ref . [jn| . 

The virtual corrections to the cross section can be written as 

a v (s,bg - A°b) = di 2ReY J [(M v )^M% (3.1) 

with t max = and t m i n = m 2 A — s and again the summatio^with bar means the same oper- 
ations as appeared in Eq. fl2.4p . M v is the renormalized amplitude for virtual corrections. 

After the renormalization procedure, a v is UV-finite. Nevertheless, it still contains 
the soft /collinear IR singularities 

d&0 (4jL + it) , (3.2) 



IR 



r(l-e) fAirfi: 



2vrr(l - 2e) 
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where 



A v 



17 

y 

47 

— + 31n- 

o s 



m. 



3 s 



(3.3) 



The soft divergences can be cancelled by adding with the soft real gluon emission cor- 
rections, and the remaining collinear divergences are absorbed into the parton distribution 
functions, which will be discussed in the next subsection. 

3.2 Real gluon emission corrections 

The 0(a s ) corrections to bg — ► A°b due to real gluon emission (shown in Fig. 4) give the 
origin of IR singularities which cancel exactly the analogous singularities present in the 
0{a s ) virtual corrections mentioned in above subsection. These singularities can be either 
of soft or collinear nature and can be conveniently isolated by slicing the bg — ► A°b + g 
phase space into different regions defined by suitable cutoffs, a method which goes under 
the general name of Phase Space Slicing(PPS). 




Figure 4: Feynman diagrams for the subprocess of bg — > A°bg with a real gluon emission. 
In this paper, we calculate the cross section for the 2 — > 3 process 

HPi) + 9{P2) -» A°(k 3 ) + b(k 4 ) + g(h), (3.4) 



adopting the method named two cutoff phase space slicing method [|20|| . We define the 
invariants 

s = (pi + P2) 2 , t = (pi - k 3 ) 2 , u = (pi - A; 4 ) 2 , 
s 45 = (fc 4 + k 5 ) 2 , ti 5 = (pi - k 5 ) 2 , i 2 5 = (P2 - h) 2 , U5 = {h - h) 2 , (3.5) 

and describe this method briefly as follows. Firstly, by introducing an arbitrary small soft 
cutoff 5 S we separate the 2^3 phase space into two regions, according to whether the 
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energy of the emitted gluon is soft, i.e. E§ < 5 s Vs/2, or hard, i.e. E5 > 5 s ys/2. The 
partonic real cross section can be written as 



a R (bg 



A°bg) = a s {bg -> A°bg) + a H (bg -» A°bg), 



(3.6) 



where <r 5 is obtained by integrating over the soft region of the emitted gluon phase space. 
a s contains all the soft IR singularities. Secondly, to isolate the remaining collinear sin- 
gularities from <7 , we further decompose a into a sum of hard-collinear (HC) and hard- 
non-collinear (HC) terms by introducing another cutoff <5 C named collinear cutoff 



& H {bg ^ A bg) = & RC {bg ^ A bg) + a UC {bg ^ A%g y 



(3.7) 



The HC regions of the phase space are those where any invariant £15, £25^45 becomes smaller 
in magnitude than 5 c s, in collinear condition, while at the same time the emitted gluon 
remains hard. <r HC contains the collinear divergences. In the soft and HC region, a s and 
<t hc can be obtained by performing the phase space integration in <i-dimension analytically. 
In the HC region, a HC is finite and may be evaluated in four dimensions using standard 
Monte Carlo techniques [21]. The cross sections, a s , <r HC and <7 HC , depend on the two 
arbitrary parameters, 5 S and S c . However, in the total real gluon emission hadronic cross 
section a R , after mass factorization, the dependence on these arbitrary cutoffs cancels, as 
will be explicitly shown in Sec. 4. This constitutes an important check of our calculation. 
The soft region of the bg — ► A°b + g phase space is defined by 



< < SmVS/2 



(3.8) 



The differential cross section in the soft region is given as 



with 



A s 



A? 



-^0 



da s = da 



a, r(i 



2vrr(l - 2e) 



47T/i; 



A s A s 

-2- 4- — i- 4- A S 

2 + e +A 



17 

T' 



— In S s 
3 



31n- 



m 



m 



2 ' 



— In 2 5 S + 6 In 5 S In • 



-t 



— In <5 S In — 

3 8 



m 



A 



s — m" A 
6 I 



+ -ln 2 



m 



6 S 



m 



+ 3Li 2 [- 



-u 



m 



A 



(3.9) 



(3.10) 



In the limit where two of the partons are collinear, the three body phase space is 
greatly simplified. And in the same limit, the leading pole approximation of the matrix 
element is valid. According to whether the collinear singularities are initial or final state 
in origin, we separate <r HC into two pieces 
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(3.11) 
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a^ is the cross section arising from the case that the emitted gluon is collinear to the 
initial partons, < £15, £25 < 5 c s. And <5"j C arises from the case that the emitted gluon is 
collinear to the final parton, < £45 < S c s. 



The cross section cr^ c can be written as 



daj c = da c 



a s T(l — e) ( Ann, 
2nF(l - 2e) 



A 



b^bg 



+ A b ^ b9 



(3.12) 



where 



j^bg 



A, 



b~.bg 



C F (3/2 + 21n5 s ), 

C F [7/2 - vr 2 /3 - In 2 5 S - ln£ c (3/2 + 2 hi 5,)]. 



(3.13) 



The cross section cp c can be written as 



da 



HC 



da 



a, r(i 



2vr T(l - 2e) 



2\ e 



(~~) S c [Pbb(z,e)G b / A (x 1 /z)G g / B (x 2 ) 



dz 1 — z 

+ Pgg(z, e)G g/A {xi/ z)G b/B (x 2 ) + {x\ <-> z 2 )]— (—^)~ e dxidx 2 . (3.14) 

where G bg i AB {x) is the bare parton distribution function, Pfcb(z,e) and P gg (z,e) are the 
(i-dimensional unregulated (z < 1) splitting function related to the usual Altarelli-Parisi 
splitting kernels[22]. Pa(z,e) (i = b,g) can be written explicitly as 



P u {z,e) = P ii {z) + eP' ii {z) 



Pbb(z) = C F 
P gg {z) = 2N 



1 + z 2 

r z 1 — z 

+ 



(* = b,g), 
Pi b (z) = -C F (1 - z) 



1 - z 



z 



+ z(l-z)}, 



P'ggi*) 



0. 



(3.15) 



with N = 3 and C F = 4/3. 

In order to factorize the collinear singularity of cj| ic into the parton distribution func- 
tion, we introduce a scale dependent parton distribution function using the MS convention: 



Gi/ A {x,n f ) = G i/A (x) + (--) 



as r(l - e) / 47^2. 
2vrr(l-2e) 1 \x 2 f 



dz 



Pu(z)G i/A (x/z), 



(i = b,g). (3.16) 



By using above definition, we replace G Si b/A,s m Eq.(3.14) and the expression for the initial 
state collinear contribution at 0(a s ) order is 



da 



HC 



Oi s r(l — e) ^ 47T/x: 
27rr(l - 2e) 
^4f c (a — > a#) 



,2\ e 



<icT° 

+ £ 

+ (a?i *-> x 2 )}dx\dx 2 , 



{Gg/ A (x 1 ,Hf)G b/B (x 2 ,iJ,f) + G g/A {x 1 ,fi f )G b/B {x 2 ,fi f ) 
+ a>g)]G g/A (xi, n f )G b/B (x 2 , n f ) 

(3.17) 
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where 



A?(b->bg) = C F (21n<f s + 3/2), 

At c (g gg) = 2Nhx5 s + (UN - 2n f )/6, 

Ag c = ^fln(4). (3.18) 



And 

•1-5. 



y s dy 

G a /A,B{x,Vf) = — G a/AtB (x/y,Hf)P aa (y), (a = g,b), (3.19) 

J x y 



with 



Paa(y) = P aa H5j-^A I )-P' aa (y,), (a = g,b). (3.20) 

y fj>f 

We can observe that the sum of the soft (Eq.Q3.9D), collinear (Eg. (|3 . 12j) , (|3. 1 7|) ) , and ultra- 
violet renormalized virtual correction (Eq.Q3.2j)) terms is finite, i.e., 



Ai + a\ = 0, 

Af + A\ + A h C h9 + (6 - bg) + Afte - 55) = 0. (3.21) 

The final result for the 0(a s ) correction consists of two contributions to the cross section: 
a two-body term and a three-body term . 



2vr 



J dx 1 dx 2 d&°{G g/A (x 1 ^ f )G b/B (x 2 , f i f )[Ai + A% + A b ^ b9 + A s c (b -> 6 5 ) + A^(g - 
+ G g / A (x 1 ,fif)G b / B (x 2 ,fJ,f) + G g / A {xuHf)G h / B {x2,Hf) + {xi <-> x 2 )}. (3.22) 

And 

(j( 3 ) = y dxidx2[G g / A (xi, fif)G b / B (x 2 , Hf) + (xi <-> x 2 )]d<T (3) , (3.23) 
with the hard-non-collinear partonic cross section given by 

^ (3) = ^- / El M 3( & 5 - ^<?)| 2 dr 3 . (3-24) 

Finally, the NLO total cross section for pp(ov pp) — » + X is 

£J 7VLO = a + CT (2) +a (3) > (325) 

4. Numeric results and discussion 

In the following numerical evaluation, we present the results of the cross section for the 
pseudoscalar Higgs boson production associated with a single high-p-p bottom quark via 
subprocess bg(bg) — ► A°b(A°b) at the LHC and Tevatron. At the LHC, the 6-jet is required 
to have a transverse momentum cut prip) > 30 GeV and a rapidity cut \r](b)\ < 2.5. At 
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the Tevatron, the b tagging regions are taken to be \r)(b)\ < 2 and pr(b) > 15 GeV. The 
SM parameters are taken as: m t = 174.3 GeV, mz = 91.188 GeV , mw = 80.419 GeV and 



dEW = 1/128 [ 23 1 . The factorization scale is taken as \if = tha/A and the renormalization 
scale is taken as fj, r = niA- We use the one- loop formula for the running strong coupling 
constant a s with a s {mz) = 0.117. 

The relevant MSSM parameters in our calculation are: the parameters Mq ^ ^ and 
Aft, in squark mass matrices, the higgsino mass parameter fx, the masses of the gluino m g 
and the A Higgs boson niA, the ratio of the vacuum expectation values of the two Higgs 
doublets tan /3. The squark mass matrix is defined as 




M\ = ( < a ^ ) (4.1) 
\a q m q m\ R J 

with 

m~ L = Mq + m q + m\ cos 2(3 (i| — e q sin 2 0w), 

2 2 2 2 2 

m qR = Mry ^ + m q + m z cos 2(3e q sin 9w 
a q = A q — //{cot (3, tan 0}, (4.2) 

for {up, down} type squarks. /| and e q are the third component of the weak isospin and 
the electric charge of the quark q. The chiral states qi and qn are transformed into the 
mass eigenstates q\ and q\ : 

= &=( C ° S ^ 8in ^V (4.3) 

Then the mass eigenvalues m qi and m q2 are given by 

= WM 2 q {R i ) ] (4.4) 

For simplicity, we assume Mq = My = M^ = At = A b = m g collectively denoted by 
MsuSY- 

In our calculation, we use the CTEQ5L parton distribution functions 1 25]. The MS 
bottom quark mass m b can be evaluated by using the next-leading order formula [26]. In 
the following equations, we use fn b {Q) to denote the MS bottom quark mass. 

m b (Q) = U 5 (Q,m b )m b (m b ), for Q < m t , 

rn b (Q) = U e (Q,m t )U 5 (m t ,m b )m b (m b ), for Q > m t , (4.5) 
where fn b = m b (m b ) = 4.3 GeV. The evolution factor Uf(f = 5,6) is 

d if) = J(f) = 8982 - 504/ + 40/ 2 

33-2/' 3(33 -2/) 2 { ' 
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In the supersymmetry limit, bottom quarks only couple to the neutral Higgs doublet H®. 
However, supersymmetry is broken and the bottom quark will receive a small coupling to 
the Higgs doublet H® from radiative corrections. Considering these corrections, ^4° Higgs 
boson coupling to the bottom quarks is given as [27] 



g w tan/? 
- m b^rz — 75 



m b g w tanf3 



2m\v 1 + A fe 2m\v 

The explicit form of A& at one-loop is given by [28, 29, 3C] 



"75 



(4.7) 



as -rrig tan &I(m~ bi , m~ h2 , m- g ) + -^A t fi, tan /3I(m fl , m~ t2 , /i), 



3vr 



(4.8) 



where ^ = u 

1 u sin p ' 



/(a, 6, c) 



2AQGeV and the function 7 is given by 
1 



(4.9) 



Fig. 5 shows that our NLO-QCD result does not depend on the arbitrary cutoffs 5 S 
and S c by using the two cutoff phase space slicing method. The two-body(cr( 2 )) and three- 
body(cr( 3 )) cross sections and the NLO cross section (a NLO ) at the LHC, are shown as the 
functions of the soft cutoff S s with the collinear cutoff 5 C = 5 s /50 . The supersymmetric 
parameters are taken as \i = tua = 200 GeV, Msusy = 500 GeV and tan (3 = 4. We can 
see the NLO cross section 

a NLO ig 

independent of the cutoffs. In the following numerical 
calculations, we take S s = 10 -4 and 5 C = 5 s /50. 
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Figure 5: Dependence of the cross sections for the A°b production at the LHC on the cutoff 5 S 
with S c = S s /50. 
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Fig. 6 shows the dependence of the LO and NLO cross sections of the process pp 
(or pp) — > bg(bg) — > A°b(A°b) + X at the LHC and Tevatron on the mass of A Higgs 
boson(mA). Here we take fi = 200 GeV, M SU sy = 500 GeV and tan/5 = 4, 15, 30. The 
relative NLO-QCD corrections are about 40% ~ 50% at the LHC and about 45% ~ 80% 
at the Tevatron when itla varies from 200 GeV to 800 GeV for all the values of tan (3 we 
have taken. 




1E-5-| , 1 , 1 , 1 , 1 , 1 , 1 , 1 1 | , 1 , 1 , 1 , 1 , 1 , 1 , 1 ME-3 

100 200 300 400 500 600 700 800 100 200 300 400 500 600 700 800 

m A [GeV] m A [GeV] 



Figure 6: The dependence of the cross sections of process pp (or pp) — > bg(bg) — > A°b(A°b) + X 
on the vtiA with tan/3 = 4, 15, 30 at the LHC and Tevatron. 

Fig. 7 shows the cross sections of the process pp (or pp) — » bg(bg) — > A°b(A°b) + X 
at the LHC and Tevatron as the functions of the ratio of the expectation vacuum values 
tan/3. We take fi = 200 GeV, M SU sy = 500 GeV and m A = 200, 500, 800 GeV. 
Since the coupling between A Higgs boson and bottom quarks is greatly enhanced with 
large tan/?(see Eq.(|2.3|), Eq.([4.7D), the cross section of process pp (or pp) — > bg(bg) — > 
A°b(A°b) + X at the LHC and Tevatron can be rather large. We can see that the cross 
section of the process pp (or pp) — > bg(bg) — > A°b(A%) +X can reach dozens of pico bar(or 
about 1 pb) when tjia = 200 GeV nad tan (3 = 40(or when = 500 GeV and tan [3 = 40) 
at the LHC, while the cross section can be 200 fb when txia = 200 GeV and tan j3 = 40 at 
the Tevatron. 

In order to study the decoupling behavior of SUSY QCD correction, we push Msusy 
to a large value. The relative SUSY QCD correction is defined as 

6a SQCD 

Asqcd = — , (4.10) 

where 5a S( ^ CD is the cross section correction contributed by the SUSY QCD diagrams 
shown in Fig. 3. In Fig.8, we depict the relative SUSY QCD correction A$qcd as the 
functions of Msusy at the LHC. The relative SUSY QCD correction is small (about 2.2% ~ 
2.5%) but not vanishing with Msusy up to 2 TeV . 
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Figure 7: The dependence of the cross sections of process pp (or pp) — > bg(bg) — ► ^4°6(A°&) + X 
on tan/3 with m A = 200, 500, 800Gey at the LHC and Tcvatron. 
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Figure 8: The dependence of Asqcd of process pp — > bg(bg) — > A°b(A°b) + X on Msusy at the 
LHC. 



To analyze the scale dependence of the cross sections, We introduce the ratio of the 
cross section at scale fi and the cross section at scale [i = tua and depict the a(fi)/a(fi = 
tua) as a function of n/rriA at the LHC in Fig. 9. For the solid lines, we fix the renormaliza- 
tion scale fj, r = tua and only show factorization scale [if dependence of the cross sections. 
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For dashed lines, (if and fj, r are taken to be identical, vary from myi/5 to 2mA- The scale 
variation of the NLO-QCD cross section may be serves as an estimate of the remaining 
theoretical uncertainty of the high order corrections. Fig.9 shows that it is evident that 
the one-loop NLO-QCD corrections reduce the LO scale dependence for both the cases. 




Figure 9: The variation of the o{\i)j(j{\i — tua) with the ratio fi/mA of process pp — > bg(bg) — ► 
A°b(A°b) + X at the LHC. 

In summary, we have computed the production of pseudoscalar Higgs boson A asso- 
ciated with a single high-pr bottom quark via subprocess bg(bg) — > A°b(A°b) including the 
NLO-QCD corrections in the MSSM at the LHC and Tevatron. We find that due to the 
enhancement of the Yukawa coupling strength of the bottom quarks with ^4° Higgs bosons 
at large tan /?, the cross section of the pp (or pp) — > A°b(A°b) can reach dozens of pico 
bar at the LHC and hundreds of fermi bar at the Tevatron. The NLO-QCD corrections 
vary between 40% ~ 50% at the LHC and 45% ~ 80% at the Tevatron respectively, in the 
parameter space we have chosen. 
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